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The value of the bottom quark mass at Q = Mz in the DR scheme is an important input for 
the analysis of supersymmetric models with a large value of tan Conventionally, however, the 
running bottom quark mass extracted from experimental data is quoted in the MS scheme at the 
scale Q = mb. We describe a two loop procedure for the conversion of the bottom quark mass from 
MS to DR scheme. The Particle Data Group value m^^{ml'^^) = 4.2 ± 0.2 GeV corresponds to a 
range of 2.65-3.03 GeV for mp^(Mz). 
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I. INTRODUCTION 

Supersymmetry (SUSY) is one of the best-studied ex- 
tensions of the Standard Model (SM) Despite the 
lack of direct evidence, many theorists consider the idea 
promising enough to examine even quantum corrections 
in SUSY theories: these include corrections to SM quan- 
tities due to sparticle loops as well as corrections to 
processes involving the production and decays of super- 
symmetric particles Quantum corrections in SUSY 
theories are also important for obtaining the implications 
of top-down models, where simple assumptions about 
symmetries of physics at the high scale result in a highly 
predictive scenario. The best known example of these is 
the correction [Q to the relation between the down type 
mass and the corresponding superpotential Yukawa cou- 
pling, which is crucial to include when analysing SUSY 
models for large values of the parameter tan /3. 

For large tan/? values, third generation sparticle and 
Higgs boson masses as well as sparticle decay patterns, 
and thus collider signals for SUSY, are sensitive to the 
value of the bottom (and to a lesser extent tau) Yukawa 
coupling 1^ . Knowledge of Yukawa couplings is also nec- 
essary for an analysis of the interesting possibility that 
these are unified at a high scale. Yukawa coupling uni- 
fication may, for instance, be realized in the minimal 
50(10) SUSY model provided tan/5 is large §. The in- 
terpretation of super-Kamiokande atmospheric neutrino 
data 1^] in terms of neutrino oscillations originating in a 
non-trivial flavour structure of the neutrino mass matrix 
has led to considerable recent interest ^ in such a 
scenario. 

To include radiative corrections in the phcnomenolog- 
ical analysis of supersymmetric models, the DR regu- 
larization by dimensional reduction pO[, a modification 



of the conventional dimensional regularization is com- 
monly used as a candidate for an invariant regulariza- 
tion of supersymmetric theories. Since the bottom quark 
Yukawa coupling is determined from the running bottom 
mass parameter at low energy, a careful determination of 
mi,{Mz) in the DR scheme is necessary. 

On the other hand, it became a standard practice to 
use experimental data to extract the value of the run- 
ning bottom mass parameter mb{mb) in the MS renor- 
malization scheme. The most recent results based on 
NNLO QCD analyses of the T sum rules can be found 
in pd[ | (see also the summary in [ p2[ which, however, 
also includes older (and by now obsolete) results in its 
world average). It is useful to explicitly examine what 
the allowed range of m^^(mb) jl^ translates into when 

expressed in terms of m^^{Mz). While an immediate 
motivation to study this is because the allowed range 



of 



,DR 



[Mz) yields an estimate of the "error" within 



which Yukawa couplings in SUSY models with Yukawa 
coupling unification should be required to unify, our re- 
sult for m]^^^{Mz) would be of relevance for any SUSY 
analysis with large tan/3. 

Many authors start with the value of the bottom 

quark mass mjj'°'° to extract (A/z), either by directly 
relating the pole mass to the running mass at the high 
scale (as described in Sec. II), or by first using the pole 
mass to derive mj^^ at the scale of the bottom quark mass 
and then evolving it to Mz using renormalization group 
equations. Since it became clear in recent years that the 
pole mass of the quark suffers from infra-red sensitivity 
which precludes its accurate determination from physical 
observables, this approach leads to significant ambigui- 
ties in the resulting value of m^^{Mz). In this paper, 
we suggest a better approach to the determination of 
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tDR 



[Mz), which by-passes the use of the 5-quark pole It is possible to invert Eq. (|^), expressing 



pole 



mass. The idea is to employ recent accurate determina- 
tions of mf,(mh) in the MS scheme use the renormal- 
ization group running up to the scale ^ = Mz and use a 
relation between the MS and DR masses at that scale to 
obtain m^^{Mz)- We will show that the relation that 
connects the masses in the two schemes at the same scale 
/i exhibits excellent convergence properties thereby mak- 
ing the suggested procedure very accurate and practically 
insensitive to higher order corrections. 



II. MS AND DR MASS FORMULAS 

Our final goal will be a relation between the MS and 
DR masses. However, since in the literature the two 
masses are usually related to the pole mass, we begin 
by collecting formulae that relate the pole mass of the 
bottom quark with its running mass (at any scale) in 
both the MS and DR schemes. 

Currently, the relation between the MS and pole quark 
mass is known through order 0[al) |l|, |l7[ |2|, 
In this paper we will only use the two loop relation 
between them: 
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where as(/x) = as{p)/TT is the MS coupling constant in 
the theory with five quark flavors, L = log(m^(^)/^^), 
and A is a correction due to light quark mass effects (aris- 
ing from the two-loop diagram with a second fermion- 
loop), given by 0, 
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S{x) = —x{l + x^)-x^ + 0{x^). (3) 
6 



terms of {mb), and then use this together with 
Eq. ^ choosing = to obtain m°^(Mz). This is 
unsatisfactory for two reasons. First of all, the relation 
between and m^^{Mz) involves not so small loga- 

rithm {CpasiMz) log(M|/m^) ^ 0.3) which leads to un- 
comfortably large corrections. Second, the perturbative 



corrections to the relation between m 



pole 



and TO^^ 



[rrib) 



are large. For example, including the 0{a'l) term in 
Eq. (P leads to a correction comparable to the contribu- 
tion of the 0(o?g) term. By the same token, when invert- 
ing Eq. (|l|) to determine the value of m^°^° out of m^^, 
it makes a significant numerical difference if one expands 
the inverse equation in powers of as(/i) or just keeps the 
r.h.s. of Eq. (|^) in the denominator. These problems 
show that, as a consequence of large perturbative cor- 
rections, the final ntimericaZ results for m^^iMz) when 
obtained from Eq. (^, are not stable against higher or- 
ders perturbative corrections. Analyses that use the 
pole mass as an input for extracting the short-distance 
mass, and via this, quark Yukawa couplings, will reflect 
this ambiguity. In what follows we suggest a procedure 
wherein the problems just mentioned are considerably 
ameliorated. 



III. DR MASS FROM THE MS BOTTOM MASS 

The first of the above problems can be by-passed be- 
cause an analytic solution of the renormalization group 
equation for the running mass in the MS scheme is 
known. This was originally obtained at three loops 
and recently the four loop term has also been computed 
p^ . For the case of five active quark flavors, and up 
to NNLO renormalization group evolution, this takes the 
form 
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In the argument of A the ratio of quark masses of the 
light quarks to the bottom quark pole mass appear. For 
< X < I the exact result is approximated by the three 
terms given in Eq. (|^) within an accuracy of 1%. 

The relation between the pole mass and the quark run- 
ning masses in the DR renormalization scheme is known 
through 0{a1) [p2). The result reads: 



-^£- = l + a,(Ai)f- 



2159 



395 ^ 

Li -I- 

72 288 



— + — + — n2 
6 9 9 



(4) 



where {jji) = (/x) /tt is the strong coupling constant 
in the DR scheme and Li = ln(mP^(^)^//i^). 



where the integration constant ttims is tbe renormaliza- 
tion group invariant bottom mass. We do not need to 
know Toms because, if we denote the right hand side as 
TTT-MsFbifJ'), the running mass at scale Mz can be calcu- 
lated from a given running mass at scale mb{mb) using 
the expression. 
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Although Eq. d) has been written in the MS scheme, 
we observe that the leading term, determined entirely 
by the first coefficient of the /3-function /3o and the first 
coefficient of the quark mass anomalous dimension 70, is 
scheme- independent. At Q = Mz the scheme-dependent 
subleading terms contribute ~ 8%. Most of this comes 
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from the next to leading a^/Tr term while the three loop 
term contributes just 0.2%. The coefRcient of this two 
loop term can be readily evaluated in the DR scheme 
using, 



Po 70 
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where 71 is the coefficient of the 0{a^^) term in the quark 
mass anomalous dimension. At two loops, we obtain for 
the running bottom mass in the DR scheme. 
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To compute the bottom running mass at rUb we need 
the value of as{mb) that corresponds to the experimen- 
tal measurement, as{Mz)- To this end we will use the 
three loop analytical formula for as ||2^ in the MS 
scheme, that is the solution of the corresponding Callan- 
Symanzyk equation. In the case of five active quark 
flavours the formula reduces to, 



^ 23^ 
2 



1 - 



348 ln(i) 
529 t 



348 
529 



1 



78073 \ 
242208 ] 



(9) 



where t = ln(/i^/A^). In practice, it is convenient to 
use the value of as(Mz) to first determine A, and then 
use Eq. again for the calculation of as{mb) which is 
inserted in Eq. (||) before the iteration. The four loop 
contributions to Eq. have also been calculated |2^ 
and found to be small. It was observed that the previous 
three loop analytical solution for as (fJ.) gives a very good 
aproximation to the numerical solution of the four loop 
RGE 1^ . Finally, we note that to compute as (fJ.) in the 
DR scheme we convert Q!s(/^)ms using p7| |. 
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By eliminating the pole mass from the relations between 
the pole and running masses presented in the previous 
section, we can derive a two loop relation between the 
DR and the MS bottom masses at the same scale. ^From 
and (|^), using Eq. (10) to convert the DR coupling 



constant to the MS coupling constant, we obtain: 
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This expression also appears, but is never used, in 
Ref. H. 

We were originally motivated to consider the relation 
between MS-DR masses because it was clear that the 
leading logarithm corrections are absent because leading 



contributions to the QCD beta function and quark mass 
anomalous dimension are scheme- independent. The ab- 
sence of all log ^/ mi, terms follows from the fact that 
both MS and DR are mass-independent renormalization 
schemes and that (7^^ — 7^^^)//?^^ is regular around 
as — Q | |2^ . The renormalization group equations then 
ensure that the relation between the two masses at the 
same scale can only depend on the coupling constant 
as(/i)- 



The perturbative series in Eq. (11) exhibits excellent 
convergence of the perturbative expansion. This is to be 
expected, since the MS and the DR mass are both truly 
short-distance quantities. The infrared renormalon prob- 
lem, which is the reason for large perturbative corrections 
in the relations between the pole mass of the quark and 
any of the short-distance masses, is not relevant there. 
Note that the cancellation of large corrections happens 
order by order in perturbation theory and for this reason 
one cannot improve on the obtained relation by including 
e.g. the three loop term from the relation between the 
pole and the MS mass without including a similar term 
from the relation between the pole and the DR mass. 
These considerations suggest the following two step 



procedure for extracting 



^DR, 



{Mz), starting from 



(mb). 



• First, evolve mf^{mb) io Q = Mz using (0) which 
is valid within the MS scheme. 

• Use Eq. (O) to convert m^^{Mz) thus obtained to 



^DR 



{Mz). 



In principle, we could have by-passed the first step by 
directly using the determination of m^^{Mz) from 
three jet 6-quark production at LEP, but the associated 
errors are so large that we cannot use it for our purposes. 

The value of mP^(Mz) obtained using this two-step 
procedure is shown as a function of m^^{mb) by the 
central solid line in Fig. |l|. This line is obtained, us- 
ing Eq. (0), for as{Mz)^ = 0.1172. The upper (lower) 
solid lines correspond to a choice for as that is smaller 
(larger) than its central value by 0.002. For comparison, 

we also show by the dotted line the value of m^^{Mz) 
using as(A/z)j^ = 0.1172. 

An alternative approach would be to derive mj^^ (mf,) 
using Eq. ( pi] ) and then evolve it to Q = Mz using 
Eq. (H). If we do so, we obtain the dashed lines in the 
figure. The relative difference between these and the 
solid lines is similar to the size of the three loop term 
in Eq. (||). The agreement between the dashed and solid 
lines provides additional evidence for the reliability of our 
procedure. 



IV. SUMMARY 



In summary, we have suggested a procedure to convert 
the value of m^^{mb) that is extracted from experimen- 
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FIG. 1: mb{Mz) as a function of m]^^ {mi,) . The solid lines 
show the results for m°^(Mz) using the first procedure de- 
scribed in the text, where we first run to Mz, and 
then convert to the DR mass at the same scale. The lines 
correspond from top to bottom to values of aa{Mz)^ ~ 
0.1152, 0.1172, 0.1192. The dashed lines show the correspond- 
ing results when we first convert to the DR scheme, and then 
evolve to Q = Mz- Finally, the dotted hne shows the value of 
mp(Mz) for q,(Mz)ms = 0.1172. 



tal data to m^^{Mz). This value serves as an important 
input for SUSY analyses with large tan p. 



Our main results are summarized in Fig. y. Most of 
the recent analyses indicate that m^^(mb) = 4.2 ± 

0.1 GeV. This translates to rn^iMz) = 2.83±0.11 GeV. 
If instead we take the Particle Data Group range 
4.0 - 4.4 GeV for m^^{mb), the range of mf^{Mz) ex- 
tends from 2.65 - 3.03 GeV. 
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